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"q ■ The hydrodynamics equations of binary mixtures of Bose gases, Fermi gases, 

and mixtures of Bose and Fermi gases in the presence of external potentials are 
^-j. | derived. These equations can be applied to current experiments on mixtures 

CN | of atomic gases confined in magnetic traps. 

The recent discoveries of Bose-Einstein condensation of atomic gases have generated 
much interest in studying dilute Bose gases trapped in external potentials. While most of 
the experiments are performed on trapped atoms of the same kind, the current traps can also 
hold different kinds of atoms at the same time, providing opportunities to study mixtures 
of Bose gases, mixtures of Fermi gases and mixtures of Bose and Fermi gases. In a recent 
C — ■ paper, we have shown Q that there is a great variety of ground state structures in the 
binary mixtures of Bose gases. These structures can be experimentally accessed by varying 
the number of trapped Bosons of each kind. The two Bosons can be different alkalis, different 
isotopes of the same alkali, or different hyperfine spin (F) states of the same alkali isotope. 
Recently, a mixture of different hyperfine states of 87 Rb (F = 1 and F — 2) has been 
produced by the JILA group ||, creating a system of interpenetrating Bose superffuids for 
the first time. 

In the case of superfluid 4 He, its superfluid properties are contained in the well known 
two-fluid hydrodynamic equations of Landau [§] that incorporate the effect of broken gauge 
symmetry into hydrodynamics. In order to understand the properties of mixtures, it is 
important to derive the corresponding hydrodynamic equations, which is the goal of this pa- 
per. Our derivation is a standard one. It follows from general principles such as conservation 
laws, broken symmetries, thermodynamic relations and the second law of thermodynamics 
It will be clear from the derivation that the form of the hydrodynamic equations is 
independent of the density of the g clSCS, clS long as the symmetry of the system remains 
unchanged. What distinguishes a dilute and a dense system is the explicit form of the co- 
efficients entering these equations, which are related to the thermodynamic potential of the 
system. In this paper, we will focus only on the derivation of the hydrodynamic equations 
of mixtures. The solutions of these equations require the evaluation of the thermodynamic 
potentials which will be discussed elsewhere. (We have recently solved the hydrodynamic 
equations of single component dilute Bose gases in spherical traps and have obtained all 
the first and second sound modes (classified by radial and angular quantum numbers) as a 
function of temperature M). 

Before proceeding, it is necessary to comment on the broken symmetry of the alkali 
gases in magnetic traps. Since the alkali atoms carry spin, the order parameter in the Bose 
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condensed phase should be a spinor field instead of a scalar field. However, in a magnetic 
trap, only the spin states locally aligned with the magnetic field B(r) are trapped |J. The 
resulting quantum gas therefore has the spin field "frozen" along the local field, leaving 
particle motion as the only degree of freedom. As a result, the spin field reduces to a scalar. 
On the other hand, even though the spin field is frozen, the spinor nature is not completely 
erased. In a recent paper |J, we have shown that the spatial variations of the spin field will 
give rise to a Berry phase which can lead to a number of dramatic effects. While the Berry 
phase effects in current traps are small, they can be magnified by changing the geometry of 
the trap. For simplicity, we shall only consider those traps where the Berry phase effects are 
negligible. One should, however, keep in mind that the hydrodynamics equations derived 
here are limited to traps of this kind. 

Since both scalar Bose superfluids and s-wave Fermion superfluids have identical broken 
£7(1) gauge symmetry, and since the hydrodynamics of a system depends entirely on the na- 
ture of its broken symmetry, the hydrodynamic equations of these two systems are identical. 
The hydrodynamics of binary mixtures of quantum gases (which can be either Fermions or 
Bosons) therefore fall in three categories, namely, mixtures of superfluids, mixtures consist- 
ing of a normal fluid and a superfluid, and mixtures of normal fluids. The superfluids can be 
the scalar condensates of Bosons or the s-wave BCS condensates of Fermions. The normal 
fluid can be the either Bosons or Fermions in their normal state. 

Among all the cases considered, the hydrodynamics of binary superfluid mixtures in a 
trap is the most complex. We shall focus on this case. As we shall see, once we derive the 
hydrodynamic equations for this system, extensions to other mixtures is easy. For a binary 
mixture of Bosons, in general, there will be regions in the trap where the two species coexist 
and regions where there is no coexistence. Furthermore, in the region of coexistence, each 
Boson species may be in their normal or superfluid states. Since the hydrodynamic modes 
in each of different regions are connected to one another, not only does one need to know the 
hydrodynamic equations of interpenetrating superfluids, but also those of the surrounding 
mixtures of normal and superfluids. 

The organization of this paper is as follows. In Section I, we will discuss the hydrody- 
namic variables of binary mixtures of scalar superfluids and their Galilean transformation 
properties. The hydrodynamic equations for this system of interpenetrating scalar super- 
fluids will be derived in Section II. The explicit form of the hydrodynamic equations (with 
and without dissipation) will be presented. Section III is a summary of the hydrodynamic 
equations for other mixture systems. 

I : The hydrodynamic variables of a binary superfluid mixtures, conservation 
laws, and their Galilean transformation properties : 

Before studying hydrodynamics, it is necessary to study the equilibrium thermodynamics 
of homogeneous systems, which is described by the entropy function S. S is a function of 
a set of thermodynamic variables {Aj}, which are quantities specifying the physical state 
of the system, S = S({Aj}). To study binary mixtures, it is useful to consider the entropy 
functions of the following sequence of systems with increasing complexity : 
(i) Single component normal fluids : In this case, {Aj} = (E, P,N,V), S = S(E, P, N, V), 
where E,P,N and V are the total energy, momentum, particle number and volume of the 
system respectively. The meaning of E, P and V will remain unchanged in the following 
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examples (ii) to (v). 

(ii) Binary mixtures of normal fluids : Each component in the mixture can be either a Boson 
or a Fermion. In this case S = S(E, P, N\, N 2 , V), where Ni and -/V 2 are the particle numbers 
of the two normal fluids labelled as 1 and 2 respectively. The reason that S is a function of 
Ni and A2 but not a function separate momenta Pi and P 2 is because while Boson 1 does 
not turn into Boson 2, they can exchange momenta through particle collisions. As a result, 
the equilibrium state of the system is characterised by a single momentum, 
(hi) Single component scalar superfluid : This is the well known case of 4 He. The system 
now acquires a broken gauge symmetry < ^(r) >= | < -0(r) > \e ie<yr \ where tjj is the field 
operator. In this case, S = S(E, P, N, V; v s ), where v s = {Ti/m)V6 is the superfluid velocity 
describing the phase winding in the system, and m is the mass of the Boson. 

(iv) Binary mixtures consisting of a scalar superfluid and a normal fluid : Denoting the 
superfluid and the normal fluid components as 1 and 2 respectively, we have for this case 
S = S(E, P, Ni, N 2 , V; v sl ). The reason that S is a function of total momentum P instead 
of separate momenta Pi and P 2 is the same as that for normal fluid mixtures. 

(v) Binary mixtures of superfluids : We denote the two superfluids as 1 and 2, and their 
masses by mi and m 2 respectively. From the discussion in cases (i) to (iv) above, we have 

S = S(E,-p,N 1 ,N 2 ,V;v al ,v 8 2), (1) 

with 

v SQ = — V0 a , a = 1,2. (2) 

m a 

From this point onwards, we shall consider only case (v). The hydrodynamic equations 
of cases (i) to (iv) can be obtained from the equations derived for case (v) by setting the 
appropriate superfluid densities (defined below) to zero. (See section (III)). 
The extensive property of S implies 

S(E, P, A 7 !, N 2 , V; v sl , v s2 ) = Vs(e, g, p x , p 2 , v s i, v s2 ), (3) 

where s = S/V, e = E/V, g = P/V are the entropy density, the energy density, the 
momentum density, and p Q = m a N a /V is the mass density of the species labelled a in 
the mixture, (a = 1,2). The arguments {e, g, p a , v sa } of the entropy density s will be 
referred to as the "hydrodynamic variables" of the system, {x{\. Defining temperature T, 
chemical potential (per unit mass) p a , "normal fluid velocity" v ra , and the vector b a as 
ds/de = l/T,ds/dp a = —p a /T,ds/dg = — v n /T, and ds/dv sa = —b a /T respectively, we 
can write 

1 v 2 2 b 

ds = -de - ■ dg - Tjrdpa - ^7' dv ««- ( 4 ) 

1 1 a=l 1 a=l 1 

Defining the pressure as P = T(dS/dV)E,p,N, and using the extensivity relation given by 
Eq.(^|), we have 

P = Ts - e + p, x px + p 2 p 2 + v n • g, (5) 
and the Gibbs-Duham relation 
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dP = sdT + P«d/i Q + g • dv„ - ^ b Q • dv SCf , (6) 



a=l 



which shows the pressure is a function of the intensive variables (T, p, a , v n , v sa ). In analogy 
to binary mixtures of classical fluids, it is useful to introduce the concentration variable 

c = Pi/p, l-c = p 2 /p, (7) 

where p = p\ + pi- If we define 

fi = // x c + // 2 (l - c), (8) 

then Eq.(|4]) can be written as 

Tds = de - v n ■ dg - p,dp - {ji x ~ P2)pdc - ^ b a ■ dv SQ . (9) 

a 

Likewise, Eq.(^|) can be written as 

P = Ts-e + pp + v n -g, (10) 
while the Gibbs-Duham relation can be cast in the form 

5P = s5T + g • 5w n + p5p - (pi - p 2 )pSc - ^ b a • 5v sa . (11) 

a 

The Galilean transformation properties of the hydrodynamic variables : As we shall see, the 
derivation of the hydrodynamic equations can be simplified by using the Galilean transfor- 
mation (GT) properties of the hydrodynamic variables Xj and their derivatives dS/dXi. A 
Galilean transformation is a space-time transformation (r — > r' = r + ut, t — > t' = t), where 
u is the relative velocity of the two inertial frames. Under this transformation, a scalar 
function F of r and t will change to F', which is a function of r' and t'. In particular, we 
have 

p' a = p a , a = 1,2 (12) 
g' = g + pu, where p = p 1 +p 2 , (13) 
e' = e + g • u + pu 2 /2, (14) 



and 



v^ a = v SQ + u. (15) 



While the validity of Eqs.(^) and (0) is easy to see, Eq.(|T^) follows from the Galilean 
transformation properties of the field operator of the species a, 
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^0=4(M)e X p(^ + ^). (16) 
The phase 9 a of < i[) a (r,t) > transforms as 

ffa = e a + r^p + a^l, (it) 

hence Eq. (|T3|) , which says that v s indeed transforms like a velocity field. Finally, we also 
note that the momentum density of the species a, written as g a , transforms as 

g Q -> g' a = ga + p Q u, and g = gi + g 2 . (18) 
Since the entropy involves counting of states, it is a Galilean invariant; we then have 

s' f (e + g • u + \p u2 )i (g + P u )> Pi, P2, (v,j + u), (v s2 + u)^ = s (e, g, p 1 , p 2 , v sl , v s2 ) . (19) 



Differentiating both sides of Eq. fljjp with respect to e and v SQ , we find that both T and h a 
are Galilean invariant, 

V = T, b' a = b a . (20) 
Differentiating both sides of Eq . (|19"D with respect to g and p a , we get 

v n = v n + u, and //' a = fi a - v n • u - ^M 2 , a = 1, 2. (21) 

This shows that v n also transforms like a velocity field. Finally, differentiating Eq. (pi]) with 
respect to u at u = gives 

bi + b 2 = g - pv n = g°, (22) 

where the superscript refers to quantities evaluated in the v n = frame. From Eq. (f2~l|) , 
we then have 

^« = Ma~^, a = 1,2, (23) 

and from Eqs.flTB]), (p~4[) , and (|21|), it is easy to see that the pressure P is a Galilean invariant 
and can also be written as 

P = Ts-e° + p° lPl + ij%p2. (24) 

Rotational symmetry : Under an infinitesimal rotation a, a vector such as g changes to 
g + a x g. Since entropy of the system is invariant under such a rotation, we have 

s (e, (g + a x g), p a , (v s i + a x v si ), (v s2 + a x v s2 )) = s(e, g, p«, v s i, v s2 ). (25) 

Differentiating both sides with respect to a, we have 

v n x g + bi x v sl + b 2 x v s2 = 0, (26) 



which implies bi x + b 2 x v° 2 = 0. The general form of b a satisfying Eq.fl26|) can be 
written as 

bj = Psl v° sl + p« - v° 2 ), b 2 = p s2 v° 2 + p(v° s2 - v°0, (27) 
where p s i, p S 2, and P are parameters which will be referred to as superfluid densities. With 



this parameterization, using Eq. fl22"|) , we have 

g° = p s iv° 1+ p s2 v s ° 2 , (28) 

g = PnYn + PslV s l + p s2 V s2 , Pn = P ~ Psl ~ Ps2, (29) 

where p n will be refered to as the normal fluid density. For later discussions, it is convenient 
to define the "diffusion" current 3 D as 

g° = bx + 3 D , g° = b 2 - 3 d . (30) 

Using the above definition of 3 D , g a can be written as 

gi = PnlV„ + p s iv si + p(v s i - v s2 ) + 3 d , (31) 

g2 = Pn2V n + p s2 V s2 + P(v s2 - V si ) - 3 D , (32) 

where the normal density of the species a , p na is defined as 

Pna = Pa — Psa, Pn = Pnl + Pn2- (33) 

From Eqs.( |31]) and (0), it is clear that the p term describes the backflow created by one 
super fluid component on the other in the mixture. 

Continuity equations and symmetry relations : The time evolution of the hydrodynamic 
variables {x{\ are governed by continuity equations and conservation laws, which are 

d t pa = -V • g«, a = 1, 2, (34) 
dtgi = -PiVi0i - p 2 Vi0 2 - VjUij, (35) 
d t e = -V • Q - gi • V0i - g 2 • V0 2 , (36) 

5tv sa = -Vx a , a = 1,2, (37) 

where Xa = ——dtO a , and <pi and 2 are the external potential experienced by Bosons 1 and 
2 respectively, Eqs.(^) to ([5fl) are continuity equations for density, momentum and energy 
respectively. Ily is the stress tensor and Q is the energy current. An argument of Martin 



et al, shows that n i7 - can always be taken to be a symmetric tensor [10]. The quantities 
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{Q, Hij, g a , Xa} = {Ji\i referred to as "currents", control the flow of the hydrodynamic 
variables {xj}. 

To complete the hydrodynamic description, we need to express the currents {J7i} in terms 
of the hydrodynamic variables {x{\ so as to form a closed set of equations. As explained in 
the next section, this can be done by studying the entropy production of the system. The 
derivation can be simplified considerably by using the Galilean transformation properties of 
{Ji\- By noting that the partial derivatives dt and V transform as d% — > df = dt — u ■ V, 
and V — > V = V under a Galilean transformation, we have 

Q'i = Qi + RijUj + Ute' + -(g u + g 2i )u 2 (38) 

= Uij + Ui(gi + g 2 )j + (gi + g2)iUj + puiUj (39) 

u 2 

Xa = Xa + U- V sa + —. (40) 

For example, Eq.fllQ]) is obtained by noting x' a = —^(d t — u-V)(9 a + m a u-r/h + m a u 2 t/2h) 
= Xa -u 2 /2 + u- v SQ + u 2 . 

Using these properties, we can write Eqs. fl38]) to (EOT) as 



Qi = Q°i+ Tl%v nj + v ni e + ~{gl + g%v 2 n , (41) 

U i:j = n°- + Vniigl + g% + (gl + g%v nj + pv ni v nj} (42) 

Xa = X°+v n -v° a + -v^. (43) 

For later use, we note that 

(V-g)° = V-g° + pV-v n , (44) 

(V • Q)° = V • Q° + U^djVm + e°(V • v n ), (45) 

(V Xa )° = V X ° + {v Q JiVv ni . (46) 

II. Derivation of the hydrodynamic equations for binary mixtures of scalar su- 
perfluids : 

According to the second law of thermodynamics, the entropy production of the system 
must be positive, 4 J d?xs > 0. Note that 

(47) 



7 



Evaluating Eq.(f47|) in the v„ = frame, we have 

T(d t s)° = - (V ■ Q° + U tJ d iVnj + e°V ■ 

+ E [-8° • + fJ. aPa V • v n + /i° V • g° + b Q • V X ° 

a=l,2 

+«)i(b a - V)v ni ]. 
Using Eq.(||), and Eqs.flU) to (||), Eq.flHD can be written as 
T(d t s)° = -V • [Q° + E (g°0 a - b aX °)] + 



(48) 



a=l,2 



+ E K^a + a£) v • g° - x° v ■ b c 



diV n j 



(49) 



Q = l,2 



By noting that d t s + V ■ (sv n )= ((9ts) + sV ■ v n is Galilean invariant and using Eq.([ 
Eq.(f4"9"|) can be written as 



T[d t s + V • (sv n )} = — V ■ Q D — U.piV nj - E Xa V • K + z/V • J 



1) 



(50) 



a=l,2 



where J D is defined in Eq . (|30"D . We introduce IT , Q D , x D , J D , and v D such that 

n° = V+ E UOi + ng. (si) 



a=l,2 



Q°= E ( b «x° - 0«g°) + Q 



a=l,2 



(52) 



(53) 



(54) 



The determination of the "currents" {Ji} = {Q, 11^-, g a , Xa} now reduces to the determina- 
tion of their "dissipative" components {jf} = {Q D , ng, J D , Xa }■ Introducting i?i and i?2 
defined below, Eq.([|7]) becomes 



d t s + V- Uv n + ^ —J =R = R 1 + R 2 , 



(55) 



Ri = -- [xf v ■ bx + X 2 v • b 2 + ng^ 



T 



(56) 



R 2 = Q D ■ V 



T 



J D -V 



-(Q d -zaJ^ 



/VT 
V ft 



VV 
~T~ 



(57) 
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The right hand side of Eq. ([55]) expresses the entropy generated at r by the "external forces" 
{fi} = { — VT, diV n j, — Vf, V • b Q }. The dissipative "currents" {Jf} are responses to the 
disturbances {f\} as the system tries to relax to equilibrium. In the linear response limit, 
{<Ji D } is linear in {fi}, and R becomes a quadratic form in the forces, which must be 
positive in order to satisfy the second law of thermodynamics. Once the relations between 
the dissipative currents {Jf} and the driving forces {fi} are determined, the hydrodynamic 
equations, given by Eqs.flSIp to ( |3"TD will form a closed set. Using Eqs.([JT]) and ( |3"2"| ) for g Q , 
Eqs.(|4|) and @ for r%, Eqs.(f43|) and fl53|) for X a, and Eqs.© and @ for Q, we can 
write Eqs.flSD to ((57]) as 

I. d tPa = -V-g a , a = 1,2, (58) 

gi = Pniv n + p s iv s i + p (v si - v s2 ) + J D , (59) 

g2 = Pn2V„, + p s2 V s2 + P (V s2 - V si ) - 3 D , (60) 

ii. d t ( gl + g2 ), = -piV^! - p 2 v,0 2 - v^nr. + ng), (6i) 

11^ = 5ijP + (pnl + Pn2) V ni V nj + Ea=l,2 Psa (tw), (v S a)j 

+ p (v s i ~ v s2 ) i (v sl - v s2 ) j , (62) 
III. d t \ sa = -V U a + p a + v„ • v sa + x D ) , a = 1,2, (63) 



IV. d t s + V ■ (*v B ) + V • - —J = i2, (64) 

where we have replaced the energy equation Eq . (|36|) by the entropy equation Eq.(|4]). For 
hydrodynamic processes with weak dissipation, the dissipative currents J D , U D , x D , Q D can 
be ignored and the resulting equations I to IV above will be referred to as "dissipation-free" 
hydrodynamic equations. For dissipative processes, however, it is necessary to determine 
the explicit form of the dissipative currents. 



Determination of dissipative currents : Note that the forces in Ri and i? 2 (eq.(^6|) and 
Eq.flBTD ) have different parity and time reversal symmetries. As a result, the currents 



(Xa, ng) and (Q D , J D ) are only linear combinations of (V ■ b a , diV n j) and (VT, V^) respec- 



tively. Since ng is symmetric, it can be decomposed as 



ng = Aij + rSij, TrA = 0, (65) 



which allows R± to be expressed as 
1 



Rl = —rr, 
T 



1 / 2 , 

-Aij f ViV nj + VjV ni - -SijV ■ v n ) + rV • v n + J2 Xa(V ■ b Q ) 



a=l,2 



(66) 



Linear response implies 



A n = -// | ViV nj + VjV ni - jj(%V • v„) . (67) 
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Coo Coi Co2 
Coi Cll Cl2 
Co2 C21 C22 




(68) 



where rj will be referred to as the "first" viscosity, and dj as the "second" viscosity matrix. 
The fact that Qj is symmetric is a consequence of the general Onsager relations for kinetic 
coefficients. More explicitly, we have 



n 



-V [ ViV nj + VjV ni - -5 i:j V • v„ ) - C00V ■ v r , 



C01V ■ bi - C02V ■ b 2 , 



and 



<0aV • V n - Ca^V ■ b 



(69) 



(70) 



=1,2 



Recall from Eq.fl27|) that bi 



PslV° sl 



P~( v °i - v ° 2 ) an d b 2 



Ps2V° 2 



P(v s2-0- Clearly, 



in order to have i?i > 0, we need to have r\ > 0, and det|£| > 0. 

To find the dissipative currents Q D and J D , we note that the expression for R2 is identical 
to the entropy production of a classical binary mixture ||11||. The expressions of these currents 



are therefore identical to those of classical fluids and are of the form 



(71) 



- vJ D = -PT* (^) - 7 T 2 (^) , (72) 

where we have again used the Onsager relation to identify the off-diagonal kinetic coefficients 
in Eqs.tfFID and flTD. Eq.fl57|) then becomes R 2 = 7(VT) 2 /T 2 +a(Vz/) 2 /T+2/5 (VT • Vv) /T, 
or, 

where k = 7 — [3 2 T /a. The positivity of i? 2 implies k > 0, and a > 0. Eqs.fl69|), (|70|), ([71]) , 
and ( |72D completely determine the dissipative currents and complete the full dissipative 
hydrodynamic equations, given by Eqs.(^) to ( |51D together with Eqs. (|69]) to and Eq. ([72]) . 

While we shall not be solving the hydrodynamic equations here, we remind the readers 
that the sound modes of the system are contained in the so-called linearized hydrodynamic 
equations, which are obtained by linearizing Eqs. (|58| ) to ( j64|) about the equilibrium configu- 
ration. In the single component superfluid case, Landau M has pointed out that it is useful 
to express all the hydrodynamic variables in terms of externally controllable quantities such 
as temperature T and pressure P. (Also see reference ||). In the mixture case, there is 
another controllable quantity, which is the concentration c = pi/p defined earlier in Section 
I. Expressing 3 D and Q D in terms of T, P and c as in Landau and Lifshitz [11], we have 



^ = - pD (vc+^VT + ^Vp), (74) 
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Q 



D 



vr 



dc 



-T 



P,T 



du s 
df 



v 



kVT, 



(75) 



P,c 



where k, D, k?, hp are the thermal conductivity, the diffusion coefficient, the coefficient of 
thermal diffusion, and the coefficient of baro-diffusion respectively; they can be written as 



K = 7 - p 2 T/a, D 



P \dc / 



T,P 



(76) 



Using the relations above, the full expression for entropy production is 



(77) 



77 / 2 

R = 2T V ^ iVnj + ^ jVni ~ 3 ijV ' V " 

+%(V-v„,) 2 + 2 £ % (V • v n ) (V • b a ) 

1 a=l,2 1 

+ E ^(V-b a ) 2 + +2 2 %(V-b Q )(V-b a 



«=1,2 



Q!,<T=1,2 



(VT) 2 Dp^z/ / k t k p 



(78) 
(79) 



III. Hydrodynamic equations of other mixture systems : 

(a) Binary mixture of a normal fluid and a superfluid : Let component 1 be the superfluid 
and component 2 be the normal fluid. Since 2 does not have broken gauge symmetry, v s2 is 
absent in the hydrodynamic equations. Rotational symmetry (Eq. (|20|)) yeilds bi x v°j = 0, 
which implies b x = p s iv sl , or p 2 = P = 0. Following exactly the same steps as before, it is 
easy to see that the the hydrodynamic equations of this system are still given by Eqs.([58|) 
to (|6"1D with the quantities (p S 2,P, v s2 ) all set to zero. Moreover, the dissipative part of the 
stress tensor U"^ as well as x? are given by Eqs.(|69D and (|70|) with ^ 2 and b 2 set to zero, 
while the expressions of the diffusion current J D and the heat current Q D are still given by 
Eqs.(|74]) and fl75|). 



(b) Binary mixture of normal fluids : 

In this case, both v s i and v s2 are absent from the equations. The hydrodynamic equations 
are given by Eqs. (|58|) to Eq. (|64]) with the quantities {p s x, p S 2, P, v s i, v s2 ) set to zero. The 
dissipative component U D is given by Eq. (|69|) with all ( a /3 = except for ( 00 , and x° — 0- 
The diffusion current and heat current are still given by Eqs.(|73[) and ([T5|). 

Acknowledgement : This work is support in part by NSF Grant DMR-9705295. 
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